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I. INTRODUCTION 

Theoretical investigations of the QCD phase diagram 
at low temperatures T <C T c and high baryochemical po¬ 
tential /ts ~ 1 GeV have recently received new impulses 
when color superconducting quark matter phases with 
large diquark pairing gaps of about 100 MeV were sug¬ 
gested within nonperturbative low-energy QCD models 
1, 3,y]. Prior to these works, color superconductivity in 
quark matter has been discussed using perturbative one- 
gluon exchange interactions resulting in small diquark 
pairing gaps of about 1 MeV only 0, |5j . More system¬ 
atic studies of the possible diquark pairing patterns have 
revealed a very rich phase structure of cold quark mat¬ 
ter, see e.g. _(§, ,'J, $>!§, Sflj. The most prominent color 
superconducting phases with large diquark pairing gaps 
are the two-flavor scalar diquark condensate (2SC) and 
the color-flavor locking (CFL) condensate. The latter 
requires approximate SU(3) flavor symmetry and occurs 
therefore only at rather lar ge q uark chemical potentials, 
H q > 410 — 500 MeV HEME!. 

The most favorable places in nature where color super¬ 
conducting states of matter are expected to occur are the 
interiors of compact stars, with temperatures well below 
1 MeV and central densities exceeding the nuclear satura¬ 
tion density by about one order of magnitude. However, 
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a closer examination of the compact star constraints on 
phases of dense matter reveals that the requirement of 
charge neutrality 0 may inhibit flavor asymmetric pair 
condensation in the 2SC phase and even influence the 
CFL phase border. As the quark chemical potential cor¬ 
responding to the highest possible central densities may 
barely reach the above constraint for the occurrence of 
strange quark matter 0: it is expected that the vol¬ 
ume fraction of a strange quark matter phase in a com¬ 
pact star cannot be sufficient to entail observable con¬ 
sequences. The two-flavor quark matter phases, on the 
other hand, may occupy a large volume [l5|, [lg, but due 
to the flavor asymmetry the occurrence of the 2SC phase 
is rather model dependent 1100 - 

Most sensitive observable consequences of quark mat¬ 
ter phases in compact stars come from the cooling evo¬ 
lution and it has been demonstrated that the occurrence 
of unpaired quarks in the core of hybrid stars leads to 
rapid cooling via the direct Urea process J3 0, 0|- 
To describe existing cooling data for isolated compact 
stars the occurrence of quark pairing in a 2SC phase is 
not sufficient and a residual weak pairing channel with a 
small gap Ay ~ 10 — 100 keV has been introduced phe¬ 
nomenologically (El- No microscopic calculation for the 
justification of this 2SC+X phase could be presented up 
to now. 

Given this fact and that the 2SC phase is very fragile 
with respect to the flavor asymmetry, the question arises 
whether there exist other pairing patterns in two-flavor 
quark matter which survive compact star constraints and 
fulfill the requirements from the cooling phenomenology 
that all quarks shall be paired with the smallest pairing 
gap being of the order of several keV. Single flavor spin- 
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1 pairs are good candidates since they are inert against 
large splittings in the quark Fermi levels for different fla¬ 
vors occurring due to electroneutrality They have been 
introduced first in Refs. |22l.l23| and their properties have 
been investigated later more in detail, see 0 and Refs, 
therein. Recently, a specific spin-1 pairing pattern - the 
so-called A phase - has been employed to sug gest a new 
mechanism for the origin of pulsar kicks |24j |. The idea 
is based on the fact that the pairing in the A phase is 
anisotropic, giving rise to a directed neutrino emission 
during the cooling process. On the other hand, there are 
still ungapped modes in the A phase, which would again 
lead to difficulties with cooling phenomenology. 

Therefore, in the present work we want to argue in fa¬ 
vor of spin-1 pairing in the color spin locked (CSL) phase 
as a possible candidate. This phase has recently been 
found to be the most stable spin -1 phase, at least at very 
high densities 0. As we will show, as long as the quarks 
do not become exactly massless, all modes are gapped 
in the CSL phase. To study this more quantitatively, 
we present a microscopic calculation within an NJL-type 
model, thus providing results for future phenomenologi¬ 
cal applications, like compact star cooling pi ], 

II. NJL MODEL FOR THE CSL PHASE 

The CSL pairing pattern we want to study is defined 
by the equality of three diquark condensates, 

(Of C 7 3 A 2 q f ) = ( qj C'y 1 X 7 q f ) = (qj C~/ 2 \ 5 q f ) = iff , 

(!) 

where qf is a quark field of flavor / with three color 
and, of course, two spin degrees of freedom. The opera¬ 
tors A a are antisymmetric Gell-Mann matrices acting in 
color space, while 7 ^ are Dirac matrices, and C = 
is the matrix of charge conjugation. Hence, each of the 
three diquark condensates belongs to the antisymmet¬ 
ric antitriplet in color space and, at the same time, rep¬ 
resents a vector component in spin-space. This means, 
both, color 5/7(3) and rotational 0(3), are broken for 
rjf 7 ^ 0. There remains, however, a residual invariance 
under a common transformation in which color and spin 
rotations are locked to each other. 

Together with the diquark condensates, we consider 
the presence of the quark-antiquark condensates 

<Lf = (9/9/> > ( 2 ) 

which are responsible for dynamical chiral symmetry 
breaking in vacuum. In the density regime we are in¬ 
terested in, these condensates are relatively small. How¬ 
ever, as we will see, it is important to keep them in the 
calculations. 

In order to investigate the properties of this condensa¬ 
tion pattern, we consider a two-flavor system of quarks, 
q = (u,d) T , and employ an NJL-type Lagrangian 


and 4-point interactions in the channels corresponding to 
Eqs. Q and 0: 
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= G ^2{dT a q) 2 + ••• , ( 5 ) 

a—0 

C qq = -H v (q^C T S \Aq T ) {q T C'y l _ l T S \Aq) + ■■■(6 ) 

Here the matrices t q , a = 1,2,3, are the usual Pauli 
matrices in flavor space while To is a unit matrix. As 
before, Gell-Mann matrices in color space are symbolized 
by A. The subscripts A and S denote antisymmetric and 
symmetric matrices, respectively, i.e., rs £ {to, ti, T 3 } 
and Aa £ {A 2 ,A 5 ,A 7 }. As usual, repeated indices are 
summed over. 

C qq and C qq should be interpreted as effective interac¬ 
tions to be used in mean-field (Hartree) approximation. 
In general, there could be further interaction terms, indi¬ 
cated by the ellipsis in Eqs. © and 0 . For instance, C qq 
could contain a scalar term in the flavor singlet color an¬ 
titriplet channel, corresponding to standard 2SC pairing. 
However, as outlined in the Introduction, we assume that 
under compact star conditions the mismatch between up 
and down quark Fermi momenta is too large to allow for 
cross-flavor condensation. In this case, the correspond¬ 
ing interaction terms do not contribute to the mean-field 
thermodynamics. In the same way, we assume that there 
is no significant condensation in any other channel except 
? 7 / and erf. 

As a specific example, we could think of C qq and C qq to 
arise via a Fierz transformation from a color-current in¬ 
teraction Ci n t = ~g{q^\ a q){q r )^ a q): corresponding 
to the quantum numbers of a single gluon exchange. In 
this case the ratio of the two coupling constants is given 
by G : H v = 1 : |. Later, we will employ this relation in 
order to constrain the number of parameters in the nu¬ 
merical calculations. At the present stage, however, we 
can leave G and H v as arbitrary constants. 

In order to derive the mean-field thermodynamic po¬ 
tential for an ensemble of up and down quarks at temper¬ 
ature T and chemical potentials g u and g d , respectively, 
we linearize C e g in the presence of the condensates 77 / 
and Of. One finds that the different flavors decouple, 
i.e., the thermodynamic potential is given by the sum 

Cl q (T, {/i/}) = f l u (T, Hu) + Sld{T, p d ) . (7) 

Defining the constituent quark masses and the diquark 
gaps 

Mf = rrif — 4:Gof , Af = AH v r]f , ( 8 ) 

and using the Nambu-Gorkov bispinors 

^ = 73 ( Cf ) ’ ^ = 73 ( 9, 9 T C') (9) 

the contributions for each flavor can be written as 


C e « = Cq + C q 


Ca 


with a free part 


(3) 

%(?>/) = 1 
n J 

r d 3 p 1 
(2t r) 3 2 

(4) 

+ T(«,. 

-m) 2 + ; 


877, 


■|A/| 2 , 


C 0 = q(i ft - m) q , 


( 10 ) 
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where the sum is over fermionic Matsubara frequencies 
u n = (2 n + 1 ) 7 tT. The inverse of fermion propagator is 
given by 


S/O) = 


i> + p / 7 ° — Mf 

— At 


A f 

$ — Hf'y 0 — Mf 


( 11 ) 


defining a set of gap equations for A f and Mf. Among 
the solutions, the stable one is the solution which corre¬ 
sponds to the absolute minimum of Qf. 


III. DISPERSION LAWS 


in Nambu-Gorkov space, where A f = Af("/ 3 \ 2 + 7 1 A 7 + 
7 2 A 5 ). Taking into account its color and Dirac structure, 
this is a 24 x 24 matrix. Explicit evaluation of the trace 
yields 


Tr In 



where 


E$l +in ^M +ln 2M\ 


T 2 


2 1 2 


T s y> 


Fj{y) = Po — (e 2 f + p 2 + |A/1 2 t 2 \/VH + l A /l 2 P 2 ) 

(13) 

and 

Gf(p) = (po — (£f ~ M/) 2 ) (po ~ ( e / + A 4 /) 2 ) 

- 10 |A/| 2 A (pg - (e/ - p/) 2 ) (pi - (s f + /i/) 2 ) 

- 33|A/| 4 B + B_ - 40|A/| 6 A + 16|A/| 8 . (14) 

Here, we have introduced the abbreviations e 2 = p 2 +M‘j 
and 

A = p 2 0 -^p 2 - p)-M 2 f , (15) 

B± = p 2 0 -^p 2 -p 2 f -M] 

± — i/l6p 4 + 297/x 2 p 2 + 561 . (16) 

OO ' 

Finally, we have to turn out the Matsubara sum. To that 
end we determine the zeros E"j. k of the polynomials Fj 
and G / with respect to p g , 

F^ipotP) = Po^Ff. 12 , (17) 
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Gf(po,p) = Y[(pI ~ E}. k ) , (18) 

k —3 

and apply the formula 

T E 111 + X l)) = X » + 2Tl14 i 1 + ( 19 ) 

n A ' 


to get the result 

n f (T,p f ) = ±(M f -m) 2 + -^-\A f \ 2 

~pj ^ 3 (% + 2 Tln(l + e-^)). 

( 20 ) 


This expression is the basis for the further analysis which 
will be performed on the mean-field level, i.e. at the 
stationary points 


6Qf 

dA 7 


= 0 , 


SQ r 

—- = 0 

SM f 


The six functions Ef^, defined in Eqs. and m 
correspond to the dispersion laws of the quasiparticle 
states of flavor /. In the above discussions we have as¬ 
sumed that we know these functions. This is obviously 
the case for Ef ; l and Ef. 2 which are directly read off from 
Eq. m - 


E f, 1,2 = \j £ } + P) + l A /l 2 T + \ A f\ 2 P 2 ■ ( 22 ) 

Here we have chosen the positive square-root, corre¬ 
sponding to the quasiparticle energy, while the negative 
square-root corresponds to a quasihole state. 

It is advantageous to define effective quantities, 

£ /,eff = P “ + ^f,eS > 

Mf,eS = Mf , 

M/.eff = p} + l A /| 2 , 

and 

M f 

A fieB = —L\A f \ (26) 

l^f,e ff 

in order to bring Eq. lt22ll into the standard form 


(23) 

(24) 

(25) 


Ef, 1,2 — \J (£/,eff T Pf,es)~ + A^ eff . (27) 

Note that, in contrast to the gap parameter Af, the ef¬ 
fective gap A f ie ff corresponds to the true gap in the spec¬ 
trum of the first quasiparticle mode. In fact, as we will 
see below, it is the smallest gap of all quasiparticle modes 
in our model. It is therefore the relevant quantity for 
cooling and other transport properties. 

Compared with Af, the effective gap A/ je ff is sup¬ 
pressed by an additional factor of the order Mf/pf. This 
implies that the mode Ef -, becomes gapless in the chiral 
limit when Mf = 0. This is the reason why the current 
quark mass m and the chiral condensate <jf should be 
retained, even if they are small. 

The four remaining dispersion relations Ef-k, k = 
3 ... 6 are the zeros of the fourth-order polynomial G / (p) 
CHI- Although, in principle, the exact solutions can be 
determined analytically, the resulting expressions are ex¬ 
tremely long and difficult to handle. In practice, it is 
therefore more convenient to perform a numerical search. 

Alternatively, we can construct approximate analytical 
solutions which have the advantage that they are more 
transparent. Anticipating that Af will be small, we ex¬ 
pand E 2 . fc in powers of | Ay | 2 and insert this into Eq. 11811 : 

6 / oo 

Gf{p 0 ,p) = n \ pI - E c /;fc 

fc=3 V j=0 



( 21 ) 


(28) 
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The coefficients Cj) k can be determined comparing the 
r.h.s. of this expression order by order with Eq. <d. 
At lowest order one recovers the free dispersion laws 
for particles and antiparticles, Cy°g 5 = (ey — /iy ) 2 and 

c /°4 6 = (£/ + M/) 2 , w hile the coefficients for the first non¬ 
trivial correction are obtained as solutions of quadratic 
equations. Therefore, the particle branch and the an- 
tiparticle branch split: 


E f-,3 ,5 - ( e f ~ M/) 2 + c /^3,5 IAH 2 + ■■ ■ > 

E j - 4,6 = ( e / + A 4 /) 2 + c /^ 4,6 l ^/| 2 + ■ ■ • > ( 29 ) 

where 


c (1) 

c /; 3.5 


c (1) 

c /; 4,6 



For illustration, we have plotted the coefficients c^) k , k = 
3... 6 , in Fig. Uas functions of \p\/fif. We have chosen a 
typical chemical potential yty = 400 MeV and a typical 
constituent mass My = 30 MeV. For comparison we also 
show the result for My = 0. Except for small momenta, 
the mass effect is rather small. 

As we will see in the next section, the gap parameter 
Ay is of the order 0.1 — 10 MeV in our model. We have 
checked numerically that Eq. < 1291 ) . neglecting terms of 
order | Ay | 4 , is sufficiently accurate in this case. We can 
make use of this to perform a closer examination of the 
excitation gaps in the spectrum: 

Since Ay is small, the smallest gaps are expected in 
the vicinity of the Fermi surface. Evaluating Eqs. PI 
and m at the Fermi momentum, one obtains 

£ 2 . 3 ! 5 (p 2 = m 2 - M 2 ) ~ ( 2 ± V2 |Ay | 2 , (31) 

where we have neglected higher orders in My/yty. (The 
antiparticle modes Ef-Afi are of course irrelevant for this 
discussion.) Hence, the excitation gaps are considerably 
larger than Ay : 6 ff, which is related to Ef. i, see Ea. Pl . 
In particular, none of the dispersion laws Ey ;3 _ 6 can be¬ 
come gapless. 

One still might worry about the fact that we have eval¬ 
uated the energy only at the Fermi momentum. Perform¬ 
ing an “exact” minimization, the position of the mini¬ 
mum becomes shifted. However, this shift in momen¬ 
tum is of the order | Ay | 2 , which leads to a correction 
in the excitation energy of the order | Ay | 4 , i.e., outside 
the range of validity of our approximation. For instance, 
for My = 0, the minimum of E 2 . 5 is located at the mo¬ 
mentum \p\ = Hf + |Ay| 2 /(2/iy). At this point one finds 
Ej. 5 = 2|Ay | 2 — |Ay| 4 /(4yt 2 ). Of course, the last term 
should not be trusted. But even if we did, we would 
find a gapless solution only if | Ay | > y/8 fif. This is not 
only outside the range of our approximations, but also 
far away from any realistic value of | Ay |. 



FIG. 1: Dispersion relation coefficients c^ 1 \ , k = 3...6, as 
functions of [p|/yty for fixed yty = 400 MeV. The thick lines 
correspond to a typical dynamical mass of My = 30 MeV, the 
thin lines to My =0. 


IV. NUMERICAL RESULTS 
A. NJL model 

As we have seen, in our model the different flavors de¬ 
couple. Therefore, the solutions of the gap equations can 
be discussed for a single flavor as functions of the corre¬ 
sponding chemical potential fif. Later, we may combine 
the solutions for up and down quarks (together with lep¬ 
tons) to obtain electrically neutral matter. In this article, 
all calculations are performed at T = 0. 

We solve the gap equations using the dispersion rela¬ 
tions PI for the modes 1 and 2 and the approximate 
dispersion relations (El for the modes 3-6, where we 
neglect the terms of order |Ay| 4 . The parameter sets 
(current quark mass m, coupling constant in the meson 
channel G, and three dimensional cut-off A) are shown 
in Tab. [I] and have been determined by fitting the pion 
mass and the pion decay constant to their empirical val¬ 
ues. They correspond to the vacuum constituent quark 
masses My = 300 MeV, 350 MeV, and 400 MeV, respec¬ 
tively. For the diquark coupling constant we employ the 
color-current relation H v = 3G/8. 


Set 

A [MeV] 

G A 2 

m[MeV] 

1 

664.3 

2.06 

5.00 

2 

613.4 

2.25 

5.42 

3 

587.9 

2.44 

5.60 


TABLE I: Parameter sets for NJL model. 


In Fig. El w ® show the dynamical quark mass and the 
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FIG. 2: Dependence of the dynamical quark mass Mf and the 
CSL diquark gap A/ on the chemical potential (if of the quark 
flavor / for different parametrizations of the NJL model. 


diquark gap as functions of (if. Obviously, with increas¬ 
ing constituent mass in vacuum, the critical chemical 
potentials /r/ jCr it for the chiral phase transition also be¬ 
come larger. For all practical purposes, yu/ >cr it also marks 
the onset of the CSL phase. 1 Here the size of the gap 
A/ (/t/,crit) increases from Set 1 to Set 3 since the dimen¬ 
sionless coupling constant H v A 2 becomes larger. 

The CSL gaps are strongly ///-dependent functions in 
the considered domain and we easily identify two regimes: 
while the asymptotic behavior for large density is quite 
similar for all the sets used (A/(/i/, max ) — 10 MeV), 
the low density region is qualitatively determined by the 
parametrization (A/(/// !Cr j t ) ~ 0.1 — 6 MeV depending on 
the set). In particular for Set 1 we obtain that the super¬ 
conducting gap has a very low value at /i/ jCr it and then 
rises by two orders of magnitude up to (if = 500 MeV 
(Fig. EJ). For the other sets the corresponding increment 
is at least one order of magnitude smaller. 

The gap functions in the low-density region are crucial 
in the derivation of equations of state for compact stars 
applications. Models with high /j./, cr it like Set 3 might 
not be able to stabilize hybrid stars configurations j25J. 

In Fig.Elwe show the results of the effective gaps A/ j6 ff 
defined in HJ for the dispersion relations m as func¬ 


1 For Set 2 and Set 3, ^/,crit is smaller than the vacuum con¬ 
stituent mass. Consequently, the density at T = 0 is zero for 
Hf < /i/ jC rit 3 and, thus, there are no quarks which could pair. 
This is not quite the case for Set 1. Here P'/,crit is larger than 
the vacuum constituent mass and, hence, there is a small inter¬ 

val below fif jC rit where the density is non-zero. In this interval, 
the attractive quark-quark interaction causes a Cooper instabil¬ 
ity leading to a non-zero CSL-condensate. However, since the 
density is very small, the condensate is tiny and practically un¬ 
detectable with our numerical methods. In any case, the exis¬ 
tence of a low-density gas of unconfined constituent quarks is an 
artifact of the model and has no physical meaning. 



FIG. 3: Effective CSL pairing gap A/ >e ff, Eq. (1261 . as function 
of the chemical potential fif for different parametrizations of 
the NJL model. 


tions of (if. As discussed in Sect. Mil these are the small¬ 
est gaps in the spectrum and thus the relevant quanti¬ 
ties for the compact star cooling phenomenology. As ex¬ 
pected, the size of A/ i6 ff is up to two orders of magnitude 
smaller than A/ (cf. Fig. El- We also observe that the 
strictly rising density dependence of the latter can par¬ 
tially be compensated by the drop of the dynamical quark 
mass Mf and the increase of (if. The only exception is 
found for Set 1, where A/, as described before, has an 
early onset and rises too strongly. Therefore, we obtain 
that A f^ e ff((if) is a function which increases by one order 
of magnitude for Set 1, is nearly constant for Set 2 and 
even decreasing for Set 3. 

We stress once again that the asymptotic behavior of 
the gaps at /// jlnax is quite similar while the behavior at 
M/.crit varies strongly for the different parameter sets. 

B. Density dependent coupling 

To obtain a reasonable description of modern cooling 
data, the pairing pattern for quark matter in hybrid star 
cores should fulfill a list of constraints, as it has been 
derived in a recent analysis (2l|. These were: 

1. All quarks need to be paired (which excludes, e.g., 
the pure 2SC phase). 

2. The smallest gaps should be in the range 10 — 
100 keV (which excludes, e.g., the CFL phase). 

3. The smallest gaps should have a decreasing den¬ 
sity dependence in the relevant domain of chemical 
potentials /i cr it < // < 500 MeV. 

We have shown in the previous subsection that the first 
two conditions are fulfilled for the CSL phase, but the 
third one is not. Therefore, to qualify the CSL phase 
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discussed in this paper as a good candidate for the phase 
structure of quark matter in compact stars from the point 
of view of cooling phenomenology, we would like to study 
here a small modification of the model which could ren¬ 
der at least the smallest of the CSL gaps a decreasing 
function of the density. 

Motivated by the logarithmic momentum dependence 
of the 1-loop running coupling constant a(p 2 ) — 

1/ ln(p 2 /AQ CD ) and by the fact that near the Fermi sur¬ 
face we have |p| ~ Pp ~ pf , we suggest to introduce a 
density dependence of the couplings according to 


C'iP’f) ^^Aqcd (a 4 /) > ^V-^Aqcd (M/)-(32) 

Here, the ratio 


r ( \ - _ ln (/V ,crit /Aqcd) 

a Q ° dUX/J a(p f , crit) ln(M//A QCD ) 


(33) 


relates the running coupling constant for pf to the value 
at pf tC rit- Of course, this formula, which has to be 
taken as very schematic, cannot be used at arbitrary 
small chemical potentials. In our calculations, we as¬ 
sume that it is valid for Pf > Pf, c rit and take G and H v 
at pf = pf t crit from Table |I] Aqcd is the QCD momen¬ 
tum scale which is not well known for the case Nf = 2. 
We will consider it as a free parameter to be varied in 
the limits Aqcd = 210 — 300 MeV. 



IT [MeV] 


FIG. 4: CSL pairing gap A f(pf) for different density depen¬ 
dent coupling functions G(pf), H v (pf) as defined in Eq. I .'-{21 
for Set 2 (pf, c rit = 350 MeV). The solid line represents the 
case for constant H v = 3G/8. The dashed and dash-dotted 
lines are for density dependent coupling for different values of 
Aqcd. 


In Fig. 0] we show the CSL gaps for different /in¬ 
dependent coupling functions for the Set 2 with p cr j t = 
350 MeV. We see that the pairing gaps for Aqcd = 
280 — 300 MeV are decreasing with pf. The correspond¬ 
ing effective gaps shown in Fig. 0 change by more than 
two orders of magnitude in the relevant range of chemical 
potentials, as it is required from cooling phenomenology 
calculations |2li |. 
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FIG. 5: Same as Fig. 2] for the effective CSL gaps 


V. CONCLUSION 


We have employed an NJL-type model to study a spin- 
1 single-flavor color superconducting phase which results 
from a color-spin locking (CSL) pairing pattern in two- 
flavor quark matter. This CSL phase is particularly inter¬ 
esting for compact star cooling applications since it may 
survive under charge neutrality constraints implying a 
sufficiently large mismatch between up- and down-quark 
chemical potentials to destroy the otherwise dominant 
scalar diquark condensate of the 2SC phase. 

In the specific model we used, the different flavors de¬ 
couple and could be studied separately. It turned out 
to be important to retain non-zero quark masses, even 
if they are small, because the size of the smallest gap 
in the spectrum is controlled by the value of the density 
dependent dynamical quark mass. 

Hybrid star cooling phenomenology has suggested ear¬ 
lier that quark matter in a compact star should be in a 
color-superconducting state with all quarks gapped. The 
smallest gap in the spectrum has to be in the range of sev¬ 
eral keV and should decrease with density. Introducing 
a density dependent coupling, motivated by the running 
of the QCD coupling constant, we could show that CSL 
effective gaps can be brought in agreement with these re¬ 
quirements. It would be interesting to see, whether this 
result can also be obtained within a more microscopic 
approach, e.g., using Dyson-Schwinger formalism. 

Of course, before applying these results to compact 
stars, we must construct electrically neutral quark mat¬ 
ter in beta equilibrium. However, since there is no pair¬ 
ing between different flavors, this does not lead to major 
complications. In particular, the magnitude of the CSL 
gaps should not change considerably. Therefore, since 
the latter nicely matches the requirements from the cool¬ 
ing phenomenology, one is encouraged to consider the 
CSL quark matter phase as a viable candidate for the 
superdense matter in the interior of compact stars. 
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Another interesting question concerns the possible 
electromagnetic Meissner effect 0, f26| in the CSL phase 
0 in response to a strong magnetic field typical for 
neutron stars. A more detailed discussion of this issue is, 
however, beyond the scope of the present work. 
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